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PRELIMINARY  STUDY  OF  BIAS 
AT  SELECTED  SOVIET  SEISMIC  STATIONS 


Alan  S.  Ryall,  Jr. 


Abstract 


Magnitude  wee  determined  for  five  eerthqnnkes  on  35  sad  37  May  1080,  from 
recordings  at  eight  Soviet  aeiemic  atations  on  a  4,800  km-iong  profile  from  eastern 
Kasakh  to  eastern  Siberia.  The  records  were  haad>digitiied,  and  magnitndes  were 
determined  from  traces  corrected  for  instrument  response  as  well  as  the  uncorrected 
traces.  Four  of  the  earthquakes  were  at  Mammoth  Lakes,  California,  in  the  western 
Great  Basin;  the  fifth  was  at  Tonga  in  the  south  Pacific.  Magnitude  residuak  with 
respect  to  network-averaged  m^’s  listed  in  the  ISC  Bulletin  were  positive  for  stations  at 
Yakutsk  and  Seymchan  in  Sib^a,  negative  for  raypaths  through  the  Baikal  rift  sone, 
and  slightly  positive  for  a  station  at  Semipalatinsk,  about  100  km  from  the  East  Kasakh 
test  site.  A  comparison  of  tabulated  magnitude  residuals  for  Soviet  seismic  stations  from 
previous  work  of  Ringdal  (1985),  North  (1970)  and  Vanek  et  uL  (1978,  1980)  shows 
excellent  agreement  between  these  studies.  Our  results  were  slightly  more  scattered  but 
in  good  agreement  with  previous  work.  Recalculation  of  for  88  events  recorded  on 
granite  at  the  Nevada  Test  Site  provided  a  determination  of  magnitude  bias 
(dm^  °  — 0.10±0.35)  for  the  NTS  granite  site  with  respect  to  ISC  magnitndes. 
Comparison  of  this  value  with  Sm^  for  the  Soviet  stations  provides  a  direct  measure  of 
the  magnitude  bias  of  the  NTS  area  relative  to  areas  in  the  USSR  in  which  the  stations 
are  located.  This  bias  reflects  only  difierences  in  attenuation,  and  does  not  account  for 
diflerencea  in  coupling,  focusing,  tectonic  release,  and  instrumental  effects.  If  focusing 
and  instrumental  effects  are  negligible,  the  magnitude  bias  due  to  attenuation  differences 
between  the  NTS  granite  site  and  the  station  at  Semipalatinsk  is  -0.30. 


Introduction 

On  25  and  27  May  1980  three  6-1-  earthquakes  occurred  in  the  western  Great 
Basin  near  Mammoth  Lakes,  California,  about  200  km  NW  of  the  Nevada  Test  Site.  As 
part  of  our  investigation  of  these  earthquakes  we  requested  seismograms  from  several 
seismic  stations  in  the  Soviet  Union,  including  a  station  at  Semipalatinsk,  approximately 
100  km  NE  of  the  eastern  Kasakh  test  site.  This  note  summarises  determinations  for 
five  events  on  the  Soviet  records  —  the  three  large  Mammoth  Lakes  events,  a  smaller  (m^ 
5.3)  Mammoth  Lakes  shock,  and  an  earthquake  that  occurred  at  Tonga  (m^  6.0)  during 
the  time  frame  of  the  recordings.  The  records  were  searched  for  eight  more  events  on  the 
ISC  list,  but  none  of  these  were  recorded  at  the  Semipalatinsk  station  and  most  were  not 
recorded  by  the  other  stations.  A  moderate  (m^  5.4)  Mammoth  Lakes  shock  at  16:49 
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GMT  on  25  May  and  an  earthquake  in  the  Kurile  Islands  (m^  4.8)  were  recorded  but  not 
analysed.  Table  1  summarises  information  on  the  events  used  in  this  study;  the  Kurile 
event  is  included  for  completeness. 

Data 

The  data  included  recordings  from  eight  seismic  stations  on  a  northeast-trending 
profile  from  Ksyl-Agach  in  the  southern  part  ^  the  Kasakh  Fold  System  to  Seymchan  in 
the  Northeast  Siberian  Fold  System.  Table  2  lists  the  station  coordinates  and  Figure  1  is 
a  polar  projection  of  Asia  showing  the  station  locations.  Table  3  gives  epicentral  dis¬ 
tances  and  asimuths  for  the  five  events. 

Figure  2  is  a  polar  plot  of  the  world  centered  on  a  point  (81.5  *  N,  162.0  *  W)  about 
halfway  between  the  Mammoth  Lakes  earthquakes  and  station  SEM,  showing  the  epi¬ 
center  and  recording  stations.  Figure  3  shows  the  porition  of  the  Soviet  stations  on  a 
lower-hemisphere,  equal-area  projection  of  the  focid  sphere  for  the  three  largest  Mammoth 
Lakes  earthquakes,  together  with  projections  of  the  fault  and  auxiliary  planes  for  these 
events  from  a  study  of  long-period  P-  and  surface-waves  (Given  et  tU.,  1982).  Based  on 
point-source  theory  (Keilis-Borok,  1950)  and  the  pontion  of  the  stations  on  the  focal 

I  Table  1.  List  of  Events  I 


Date 

Time,  GMT 

Latitude 

Longitude 

Depth 

800525 

16:33:44.7 

37.596 ‘N 

118.830  *W 

6.5 

6.1 

6.1 

800525 

19:44:51.1 

37.547  *N 

118.826  *W 

5.0 

5.6 

6.0 

800525 

20:35:48.5 

37.616 ’N 

118.847 ’W 

6.1 

5.3 

5.7 

800525 

23:22:18.0 

46.680  *N 

149.070  *E 

29 

4.8 

800527 

14:50:57.1 

37.472 ‘N 

118.807 ’W 

10.8 

5.7 

6.0 

800527 

13:01:37.9 

18.610  *S 

174.700 ’E 

55 

6.0 

4.8 

Table  2.  Station  Locations 


Code 

Name 

Coordinates 

Region 

BOD 

Bodiubo 

57’51’N,  114*11’E 

Baikal 

t? 

ELT 

El’tsovka 

53  *15%  86*16’E 

Altu-Sayan 

IRK 

Irkutsk 

52*16’N,104*  19’E 

Baikal 

>s 

KZL 

Ksyl-Agach 

45  *25%  78’45’E 

Kasakhstan 

SEI 

Seymchan 

62  *  53’N,  152’26’E 

Northeast 

SEM 

Semipalatinsk 

50*24’N,  80’15’E 

Kazakhstan 

i£i 

UST 

Ust’-Kan 

50*57’N,  84*45’E 

Altai-Sayan 

L 

YAK 

Yakutsk 

62*01’N,  129’43’E 

Yakutiya 

Figure  2.  Polar  projection  showing  location  of  M2unmoth  Lakes  earth¬ 
quakes  (solid  circle)  and  stations  used  in  study  (triangles).  Radius 
of  map  is  52  * . 


25  May  16=33  (1 ) 
25  May  19=44  (2) 
27  May  14=50  (3)n 


o  -  P  axis 
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Figiire  3.  Fault-plane  solutions  for  three  of  the  largest  Mammoth 
LEikes  earthquakes  (from  Given  e<  a/.,  1982).  Lower-hemisphere, 
equal-area  projection.  Solid  circles  --  location  of  raypath  to  Soviet 
stations  used  in  study. 


Table  S.  Distance  (A  *  )  and  Aibnuth  (O  *  )  to  the  Stations 


Table  S.  Distance  (A  *  )  and  Aibnuth  (O  *  )  to  the  Stations 

Sta 

0S251633  1 

05251944  1 

05252035  I 

05271450  1 

05252322  I 

05271301 

A 

e 

A 

e 

A 

e 

A 

e 

A 

e 

A 

e 

BOD 

7S.09 

334 

75.14 

334 

75.07 

334 

75.21 

334 

23.81 

311 

95.95 

329 

ELT 

86.9S 

345 

86.99 

345 

86.92 

345 

87.07 

345 

39.76 

303 

110.06 

321 

IRK 

82.93 

335 

82.98 

335 

82.91 

335 

83.05 

335 

29.24 

298 

99.15 

322 

KZL 

9S.86 

348 

95.91 

348 

95.84 

348 

95.99 

348 

47.30 

296 

114.54 

312 

8EI 

56.77 

327 

56.82 

327 

56.75 

327 

56.89 

327 

16.35 

6 

85.23 

346 

SEM 

90.78 

348 

90.83 

348 

90.76 

348 

90.91 

348 

44.24 

301 

113.66 

318 

UST 

89.43 

345 

89.48 

345 

89.41 

345 

89.55 

345 

41.37 

300 

110.82 

318 

YAK 

66.51 

331 

66.55 

331 

66.48 

331 

66.63 

331 

18.94 

331 

91.54 

337 

sphere,  P>waves  to  the  Soviet  stations  should  have  about  80%  of  the  maximum  radiated 
amplitude. 
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Figure  4  is  a  polar  projection  centered  on  a  point  (23.0*  N,  147.0  *  E)  about  halfway 
between  station  SEM  and  the  Tonga  earthquake.  The  Kurile  event  (Figure  5)  was  in  a 
poor  distance  range  (16-47  * )  for  this  study,  and  and  had  large  scatter  in  values;  it  will 
not  be  considered  in  the  detailed  discussions  that  follow.  For  the  Tonga  earthquake  a 
plot  similar  to  Figure  3  indicates  that  the  Soviet  stations  are  close  to  the  null  axis  on  the 
fault-plane  solution,  although  the  latter  is  not  well-constrained. 

Table  4  lists  instrument  parameters  for  the  86  recordings  that  were  supplied.  Values 
of  the  maximum  magnification  and  the  period  range  (corresponding  to 
V  =  0.9  V^)  were  written  on  the  records.  For  station  KZL  the  period  range  was  not 
supplied  and  the  range  given  in  the  table  was  taken  from  the  Soviet  earthquake  bulletin 
(Akademiya  Nauk  SSSR,  1983).  For  station  BOD  the  values  of  marked  on  the 
recordings  —  "2900"  for  the  NS  component  and  "2700"  for  the  Z  and  EW  —  appeared  to 
be  an  order  of  magnitude  too  low,  and  the  magnification  (52,000)  for  that  station  was 
also  taken  from  the  earthquake  bulletin.  Use  of  the  larger  magnification  for  BOD  is  sup¬ 
ported  by  Shishkevish  (1974),  who  lists  as  45,000-49,000  for  this  station  in  1970. 

For  this  study  only  vertical-component  recordings  from  SKM-3  seismographs  were 
used.  The  SKM-3  system  consists  of  a  seismometer  and  galvanometer  with  instrument 
constants  designed  to  produce  magnification  of  ground  motion  in  the  range  30,000- 
100,000  over  the  period  range  0.3-1.5  seconds.  Appendix  A  contains  a  description  of  the 
SKM-3  system,  equations  for  calculating  the  relative  amplitude  response  and  phase  shift 
as  a  function  of  period,  and  instrument  response  calculations  for  the  eight  stations  from 
which  recordings  were  obtained.  Figure  6  shows  relative  amplitude  response  for  the  eight 
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Re^onse  curves  for  SVK-M3  short-period  vertical  instru- 
ar  stations  used  in  this  study. 


Table  4.  Data  Received  for  25  and  27  May  1980 


Sta  Inst  Comp 


BOO 


IRK 

SKM-3 

NS 

17,500 

1.1  -  1.6 

SKM-3 

Z 

17,240 

1.1  -  1.6 

SKM-3 

EW 

17,430 

1.1  -  1.6 

1 

SKD 

NS,Z,EW 

1,200 

0.2  -  20 

KZL 

SKM-3 

NS 

40,900 

0.08  -  1.6 

SKM-3 

Z 

41,500 

0.08  -  1.6 

SKM-3 

EW 

40,600 

0.08  -  1.6 

SKM-3 

EW 

1,050 

0.08  -  1.2 

SEI 

SKM-3 

NS,Z,EW 

44,600 

0.2  -  1.2 

SKD 

NS,Z,EW 

1,050 

0.2  -  20 

SEM 

SKM-3 

30,050 

0.84  -  1.5 

SKM-3 

Z 

28,400 

0.8  -  1.5 

SK 

1,700 

0.41  -  10.9 

SK 

z 

1,100 

0.5  -  11 

UST 

NS,Z,EW 

50,000 

0.2  -  1.4 

SKM-3 

E-W 

5,000 

0.2  -  1.4 

SKD 

NS,Z,EW 

1,000 

0.2  -  22 

SKM-3 

SKM-3 

SKM-3 

SK-KPCh 

SK 

SK 

SK 


18,800 

37,600 

36,800 

140 

2,130 

680 

1,930 


0.3  -  1.3 
0.8  -  1.4 
0.8  -  1.4 
0.3  -  11 
0.4  -  11 
0.4  -  9.0 
0.4  -  11 
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stations,  determined  from  equations  given  by  Aranovich  et  aL  (1974).  The  Soviet  short- 
period  seismographs  peak  at  lower  frequency  (ce.  1  Hs)  than  systems  commonly  used  in 
the  US,  and  for  some  stations  (i.e.,  SEM  and  IRK  in  this  study)  the  system  gain  is  down 
by  about  a  factor  of  ten  at  10  Hi. 

Table  5  lists  the  time  period  covered  by  each  of  the  recordings.  Except  for  the  SKM 
recordings  at  station  KZL  and  the  long-period  recordings  at  UST  all  of  the  records  were 
12  hours  long.  The  KZL  record  has  twice  the  time  resolution  of  the  other  short-period 
stations  and  appears  to  be  changed  three  times  duly;  the  SKD  record  for  UST  for  25  May 
was  24  hours  long.  Unfortunately  the  records  of  primuy  interest  to  this  analysis  —  from 
station  SEM  —  were  changed  at  different  times  than  those  at  other  stations  and  overlap 
with  the  latter  for  only  two  six-hour  periods. 
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Table  6. 

Periods  Covered  by  Recordings 

Sta 

Instrument 

Time  on 
Date  Time 

Time  off 
Date  Time 

Time  on 
Date  Time 

Time  off 
Date  Time 

BOD 

25 

26 

0002 

28 

ELT 

25 

1205 

26 

0001 

■fwl 

28 

IRK 

SKMS,SKD 

25 

1207 

26 

0018 

27 

1207 

28 

K-A 

SKM 

25 

0904 

25 

1700 

27 

0905 

27 

SEI 

SKM3,SKD 

25 

1002 

25 

2200 

27 

1002 

27 

SEM 

SKMS.SK 

25 

1628 

26 

0401 

27 

0418 

27 

U-K 

SKM3 

25 

1208 

26 

0005 

27 

1155 

28 

0005 

SKD 

25 

0010 

26 

0005 

27 

1155 

28 

0005 

YAK 

SKMS,SK 

25 

1205 

26 

0005 

27 

1205 

28 

0001 

V, 

V. 


Record  Quality.  Copies  of  the  records  were  on  35  mm  Him  of  poor-to-good  qual¬ 
ity.  For  many  of  the  short-period  records  the  contrast  between  the  trace  and  the  back¬ 
ground  was  poor,  and  attempts  to  make  enlargements  of  the  waveforms  using  a  reader- 
printer  were  not  successful.  As  a  result  the  records  had  to  be  projected  on  a  viewing 
screen  and  traced  by  hand;  the  tracings  were  then  digitized  for  computer  analysis.  Qual¬ 
ity  of  the  KZL  records  was  especially  poor,  and  the  P-waveform  for  only  one  event  was 
traced  for  that  station. 

Timing.  Two  or  three  time  corrections  were  marked  on  each  of  the  records,  usually 
corresponding  to  the  be^nning  and  end  of  the  recording  period.  Clock  corrections 
changed  at  most  stations  by  less  than  a  second  per  day.  A  few  records  were  mislabeled  as 
to  sense  of  the  time  correction,  and  in  one  case  the  time  corrections  at  the  two  ends  of  a 
recording  period  were  interchanged  on  different  records  of  the  same  station.  The  largest 
clock  drift  was  0.5  second/hour  for  the  SEM  station,  but  this  rate  was  constant.  Tracings 
of  the  P-waves  were  made  from  enlarged  projections  of  the  film,  with  average  time  scale 
about  270  mm/minute.  Station  KZL  had  a  drum  speed  twice  that  of  the  other  stations, 
so  the  enlarged  traces  were  viewed  at  a  scale  of  540  mm/min.  For  most  of  the  events  the 
beginning  of  the  P-wave  (PKP  for  stations  at  A  ^  110  * )  was  identified  on  records  of  one 
of  the  better  stations  (SEI,  YAK,  BOD)  and  the  same  point  was  marked  and  timed  on 
traces  for  the  other  stations  by  overlaying  the  recordings. 

Traveltime  residuals  were  calculated  for  the  Mammoth  Lakes  earthquakes  using 
hypocentral  coordinates  and  origin  times  determined  by  the  University  of  Nevada  from 
recordings  of  a  dense  local  network,  and  traveltimes  calculated  from  the  Herrin  et  al. 
(1968)  tables.  Initially,  expected  arrival  times  were  calculated  using  the  PRED.arr  at  the 
Center  for  Seismic  Studies;  however,  this  program  calculates  times  from  the  Jefireys- 
Bullen  (1940)  traveltime  tables,  which  differ  by  more  than  three  seconds  from  those  of 
Herrin  et  of  over  the  distance  range  of  interest.  This  resulted  in  unreasonably  large 
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negative  residuals  for  most  of  the  Soviet  stations.  The  Herrin  et  aL  tables  produced  a 
pattern  of  residuals  that  was  more  reasonable  in  terms  of  known  crustal  structure  at  the 
stations  and  complex  rupturing  in  the  Mammoth  Lakes  sequence. 

The  Soviet  stations  had  average  delays  of  0.5-2.2  seconds  for  the  Mammoth  Lakes 
events,  relative  to  the  Herrin  tables  (Table  6).  These  delays  are  not  considered  to  be 
meaningful  for  the  present  study,  since  the  larger  earthquakes  of  this  sequence  appear  to 
have  been  multiple  events,  initiated  by  small  ruptures  that  recorded  at  re^onal  stations 
but  not  at  teleseismic  distance  ranges  (Given  et  at,  1982;  Lide  and  Ryall,  1984;  Wallace, 
1985).  On  the  other  hand,  traveltime  residuals  for  individual  stations  relative  to  the  aver¬ 
age  delays  for  all  the  stations  (Table  7)  are  fairly  consistent  for  the  Mammoth  Lakes  and 
Tonga  earthquakes.  Note  that  the  two  stations  north  of  Lake  Bmkal  (BOD,  IRK)  have 
early  arrivals  for  waves  travelling  southwest  across  the  Siberian  platform  from  the  Mam¬ 
moth  Lakes  events,  but  are  late  for  paths  from  Tonga  that  cross  the  Bmkal  graben.  Resi¬ 
duals  for  station  SEM  are  small,  and  their  average  is  almost  xero. 


Table  6.  Total  P-Wave  Arrival-Time  Residuals,  Seconds 


+0.8±0.3 


+0.2±0.5 


.1±0.6 


-|-0.2±0.4 


Magnitude  Datarmlnation 

Figures  7-12  show  the  di^tised  P-waves  for  the  Mammoth  Lakes,  Tonga  and  Kurile 
earthquakes  (latter  included  for  completeness  but  not  analysed).  All  traces  are  from 
SVK-M3  instruments,  and  are  normalised;  up  on  the  traces  corresponds  to  ground  motion 
up.  With  the  exception  of  KZL  most  of  the  recordings  had  good  signal-to-noise  ratio  for 
the  large  Mammoth  Lakes  and  Tonga  events.  For  the  first  (1633Z)  Mammoth  Lakes 
earthquake,  stations  YAK  and  BOD  appear  to  have  higher  frequency  content  than  some 
of  the  other  stations,  including  SEIM  (Figure  7),  but  this  could  be  due  in  part  to  differ¬ 
ences  in  instrument  response.  P-wave  spectra  for  this  event,  corrected  for  instrument 
response,  have  peaks  at  about  0.5  and  0.7  Hs;  for  station  BOD  the  peaks  are  about  equal 
in  amplitude,  but  for  the  other  stations  the  amplitude  at  0.7  Hs  is  less  than  half  that  at 
0.5  Hs.  For  the  first  Mammoth  Lakes  (1633Z)  event  the  character  of  the  P-wave  at  the 
various  stations  is  rimilar,  consisting  of  a  small  first  arrival  followed  6  seconds  later  by  a 
larger  phase.  The  time  interval  between  these  leases  is  too  consistent  for  the  second 
arrival  to  be  PcP,  since  for  the  distance  range  of  the  Soviet  stations  t(PcP-P)  should  vary 
from  60  seconds  to  lero.  The  second  phase  could  be  another  earthquake,  but  this  is 
unlikely  because  all  of  the  Mammoth  Lakes  events  have  a  large  second  arrival  at  the 
Soviet  stations.  A  more  likely  explanation  is  that  the  phase  includes  the  reflected  waves 
pP  and  and  sP.  Based  on  the  known  crustal  structure  in  the  Mammoth  Lakes  area,  a 
pP-P  time  of  6  seconds  would  be  consistent  with  focal  depth  of  about  16  km,  and  for  that 
depth  the  sP-P  time  would  be  about  8  seconds.  For  the  second  and  third  Mammoth 
Lakes  earthquakes  in  Table  1,  the  later  phase  follows  P  by  4.9  and  2.5  seconds,  respec¬ 
tively;  if  this  phase  was  pP  the  times  would  be  consistent  with  focal  depths  of  about  13 
and  6.5  km.  The  fourth  Mammoth  Lakes  event  hu  a  complex  P-signature,  and  timing  of 
the  second  phase  is  not  consistent  from  one  station  to  another. 

Magnitude  was  deternuned  from  computer  plots  of  the  digitised  P-waves  using  the 
standard  formula 

”»*  =  lo«io  +  ^(^)» 

T 

where  A  is  the  amplitude  of  vertical  ground  motion  in  nanometers  and  T  is  the  period 
corresponding  to  amplitude  A.  The  range  in  values  listed  in  the  BuUetin  of  the  Inter¬ 
national  Seiemolofiieal  Centre  suggests  that  some  observatories  measured  the  amplitude  of 
the  initial  P-wave  for  the  Manunoth  Lakes  events,  while  others  measured  the  amplitude 
of  the  second  phase.  The  lASPEI  Commission  on  Practice  Concerning  Amplitude  and 
Period  Measurement  recommended  in  1979  that  ’’the  P  wave  amplitude  measured  should 
be  that  of  the  maximum  trace  deflection,  usually  within  the  first  25  seconds  of  the  first 
onset  or  before  the  arrival  of  the  next  clear  phase.”  Other  workers  have  recommended 
measuring  A  and  T  within  the  first  few  cycles  (Zavadil,  1980)  or  even  within  the  first  3/4 
cycle  (the  ”b"  amplitude)  of  the  P-wave  (Eisenhauer,  1980).  Because  of  the  discrepancy 
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between  mnximuin  amplitude  in  the  P-coda  and  amplitude  of  the  P  onset,  we  determined 
nij  for  the  Mammoth  Lakes  earthquakes  from  the  maximum  amplitude  in  the  first  two 
cycles  and  from  the  maximum  amplitude  in  the  first  12  seconds  ~  in  both  cases  taken  as 
one-half  the  peak-to-peak  amplitude.  The  resulting  individual  and  averse  values  are 
listed  in  Tables  8  and  9,  together  with  m^nitudes  determined  by  the  International 
S^molo^cal  Centre  (ISC)  and  Moscow  (MOS),  and  deviations  Sm^  from  the  ISC  values. 

A  third  set  of  magnitudes  was  determined  from  measurements  of  A  and  T  on  traces 
produced  by  deconvolving  the  instrument  response  from  the  digitised  recordings.  The 
deconvolution  was  accomplished  using  a  program  written  by  W.  Peppin  to  calculate  the 
complex  transfer  function  of  a  seismograph,  given  the  period  and  damping  of  the 
seismometer  and  galvanometer,  together  with  the  maximum  magnification  of  the  system 
and  the  period  at  which  the  maximum  magnification  occurs.  The  records  were  tapered 
and  bandpass-filtered  (0.25-8.0  Hs)  before  deconvolution,  and  filtered  agun  (0. 1-5.0  Hs 
bandpass)  after  deconvolution.  Figures  13-17  show  P-waves  for  all  events  except  the 
Kurile  earthquake,  after  correction  for  instrument  response  and  filtering.  For  these 
records  was  determined  only  using  the  maximum  amplitude  in  the  large  phase  follow¬ 
ing  the  initial  P-wave.  Magnitudes  and  residuals  are  listed  in  Table  10. 

For  the  Tonga  earthquake,  values  of  were  determined  from  the  maximum  ampli¬ 
tude  in  the  P-wave,  from  traces  uncorrected  and  corrected  for  instrument  response.  Four 
of  the  stations  for  this  event  were  beyond  100  *  and  distance  corrections  were  taken  from 
a  curve  developed  by  Ringdal  (1985). 


Table  8.  Magnitudes  and  Residuals,  Uncorrrected  Data,  First  Two  Cycles 


05251033 

05251944 

05252035 

05271450 

Average 

Sta 

m^ 

Sm^ 

mj 

Sm^ 

m^ 

m^ 

fm^ 

BOD 

5.47 

-.63 

4.97 

-.63 

4.96 

-.34 

5.34 

-.36 

-.49±.16 

ELT* 

5.78 

-I-.08 

-|-.084:.0 

IRK 

5.53 

-.57 

4.98 

-.62 

4.84 

-.46 

■221 

-.40±.31 

KZL* 

5.58 

-.53 

-.63db.O 

SEl 

5.67 

-.43 

5.17 

-.43 

5.09 

-.21 

5.54 

-.16 

-.31±.14 

SEM 

4.73 

-.87 

4.89 

-.41 

5.80 

+.10 

-.39±.40 

UST 

4.08 

-.92 

4.92 

-.38 

5.69 

-.01 

-.60±.39 

YAK 

6.16 

Esa 

5.21 

• 

CO 

5.48 

-I-.18 

6.38 

-I-.68 

-f.l3±.44 

Avg 

5.60 

4.96 

5.03 

5.75 

ISC 

0.1 

5.0 

5.3 

5.7 

MOS 

6.3 

5.7 

5.5 

5.9 

m 


/-i 


gl633sei 


g2035sei  K  *  g2035ust 


TaUe  0.  Magaitades  and  Realdnala,  Uacomected  Data,  First  12  Secoads 


Average 

Sm, 


+.ltt  5.4« 


6.04  -.06 

6.48 

+.38 

5.81 

6.34 

+.24 

5.63 

6.18 

+.08 

1^ 

+.60 

6.10 

+.21  5.62 


+.03  5.48 


+.32  6.03 


+.18 


5 

5.51 

5.3 

5.5 

+.oe±.i2 


+.03±.10 


+.81±.07 


+.18±.10 


-.02±.08 


+.40±.ll 


Table  10.  Magaitadea  aad  Reaidaala,  Comreeted  Data,  Firat  12  Secoada 


05271450 


6.12  +.02  I  5.65 


6.34  +.24 


6.24  +.14 


5.77 

+.07 

5.99 

+.29 

SEI 

6.47 

+.87 

5.98 

+.38 

5.76 

+.46 

+.30 

+.38±.07 

SEM 

6.32 

+.22 

5.71 

+.11 

5.53 

+.23 

5.80 

+.10 

+.17±.07 

USX 

6.13 

+.03 

5.41 

+.11 

+.02±.07 

YAK 

6.74 

+.64 

6.22 

+.62 

6.10 

+.80 

6.35 

+.66 

+.68±.08 

Average 


+.07±.09 


+.23db.07 


+.14±. 


Table  11.  Magnitvdee  and  Aaaidiials  for  Tonga  Erent 


Uaeorrected  Data 

1  Corrected  Data 

Sta 

m^ 

6m^ 

Sm^ 

BOD 

5.58 

-.42 

5.BS 

-.17 

ELT 

6.09 

+.09 

6.22 

+.22 

IRK 

5.S8 

-.62 

5.66 

-.84 

SEI 

6.19 

+.19 

6.17 

+.17 

SEM 

6.07 

+.07 

6.06 

+.06 

UST 

5.79 

-.21 

5.78 

-.22 

YAK 

6.16 

+.16 

6.15 

+.15 

Avg 

5.89 

5.98 

ISC 

6.0 

6.0 

MOS 

5.9 

5.9 

Disevuaion 

Magnitude  has  been  determined  for  five  earthquakes  on  25  and  27  May  1980, 
from  recordings  at  eight  Soviet  seismic  stations  on  a  4,300  km>long  profile  from  eastern 
Kasakh  to  eastern  Siberia.  In  general,  nu^initudes  from  this  study  for  the  Mammoth 
Lakes  earthquakes  (Tables  9  and  10)  are  about  0.2  unit  higher  than  those  determined 
by  the  ISC,  but  for  the  Tonga  event  our  is  close  to  the  ISC  value.  Our  values  are 
about  the  same  as  magnitudes  attributed  to  Moscow  in  the  ISC  BuUetin:  For  the  data 
uncorrected  for  instrument  response  average  values  in  this  study  differ  from  the  Moscow 
m^  by  0.01±0.06  unit;  for  the  corrected  data  our  values  are  higher  than  those  of  Moscow 
by  0.07±0.06.  Magnitudes  detemuned  from  A  and  T  measured  in  the  first  two  cycles  of 
the  P-wave  for  the  Mammoth  Lakes  earthquakes  (Table  8)  averaged  0.33  and  0.50  m^ 
unit  smaller  than  the  ISC  and  Moscow  values,  respectively. 

The  largest  positive  residuals  were  for  station  YAK  (average  Sm^  +0.33  for  the 
uncorrected  data  in  Tables  9  and  11;  +0.42  for  the  corrected  data  in  Tables  10  and  11). 
Yakutsk  is  located  on  the  central  Siberian  platform  in  an  area  where,  according  to 
Potap’ev  et  oL  (1974),  the  depth  to  crystalline  basement  is  shallow,  0.5>1.0  km.  About  10 
km  west  of  Yakutsk  a  mqjor  north-south  fault  offsets  basement  rocks,  with  the  western 
side  downdropped  by  about  4  km.  P-wave  velocities  in  the  crust  are  high,  6.2-6.4  km/sec 
(Vol’vovsky,  1973),  the  crust  is  about  40  km  thick,  and  the  P,  velocity  is  about  8.0 
km/sec  (Vol’vovsky  and  Vol’vovsky,  1975;  Ryall  et  at,  1980). 

The  largest  negative  m^  residuals  were  observed  for  station  IRK  for  the  Tonga  event 
(-0.62  uncorrected,  -0.34  corrected),  for  which  the  raypath  crosses  the  southern  part  of  the 
Bukal  rift  sone  in  a  WNW  direction.  Station  BOD,  at  the  northeast  end  of  the  rift,  also 
has  negative  residuals  for  the  Tonga  event  (-0.31  uncorrected,  -0.12  corrected).  The  rift 


lone  is  clinracterised  by  complex  geolo^c  structure,  including  deep  crustal  inbomo- 
gmdties  and  velocity  anomalies  in  areas  of  recent  rifting.  Within  the  rift  sone  P>wave 
velodty  in  the  crust  increases  from  5.8  km/sec  near  the  surface  to  6.4  km/sec  at  the 
crust/mantle  interface,  crustal  thickness  is  34-36  km,  and  P,  velodty  is  7.7-7.8  km/sec 
(Pusyrev  et  oL,  1975;  Ryall  et  aL,  1980).  Pusyrev  et  oL  (1977)  attribute  the  low  P,  velo¬ 
city  to  partial  melting  of  upper  mantle  material,  and  compare  the  Baikal  region  with  the 
E^ast  African  rift  system,  the  RUne  graben  and  the  Basin  and  Range  province.  It  is 
interesting  that  IRK  and  BOD,  located  on  the  southern  edge  of  the  central  Siberian  plat¬ 
form  just  north  of  the  Baikal  rift,  both  have  positive  values  of  Sm^  (-H).03±0.10  and 
-|-0.06±0.12  uncorrected,  respectively;  -|-0.23±0.07  and  -H).07±0.09  corrected)  for  the 
Mammoth  Lakes  earthquakes,  for  which  the  P-waves  approach  the  stations  in  a  SSW 
direction  across  the  platform  and  do  not  appear  to  be  affected  by  the  structure  of  the  rift 
sons.  According  to  Pusyrev  et  eL  (1977)  the  platform  is  characterised  by  uniform  layer¬ 
ing,  with  distinct  reflecting  horisons  at  depths  of  18-21  and  23-27  km,  crustal  thickneM  of 
about  40  km,  and  "normal"  P^  velocity  of  about  8.1  km/sec. 

The  rcnduals  for  station  SEM  are  intermediate  between  IRK  and  YAK.  SEM  is 
near  the  Irtysh  Itiver,  on  the  boundary  between  the  west  Siberian  platform  to  the  north 
and  the  Kasakh  fold  system  to  the  south.  In  this  area  the  crustal  thickness  is  about  45 
km,  vdodty  is  8.2-8.4  km/sec,  and  P-velocity  in  the  crust  b  high  (VoPvovsky  and 
Vol’vovsky,  1975).  The  station  b  in  the  Zaysan  fold  belt,  which  b  relatively  simple  in 
terms  of  stratigraphy  and  structure.  Sediments,  principally  of  Carboniferous  age,  contain 
thick  limestone  deposits  and  lesser  amounts  of  interbedded  volcanics.  Folding  and  fault¬ 
ing  are  relatively  nunor  in  thb  area,  compared  with  more  extensive  folding,  faulting  and 
intrumon  in  the  Chingis-Tarbagatu  geanticlinal  sone  to  the  southwest,  in  which  the 
eastern  Kasakh  test  site  b  located  (Peyve  and  Mossakovsky,  1982).  For  the  Mammoth 
Lakes  earthquakes  Sm^  for  thb  station  b  -|-0.13±0.10  for  the  uncorrected  data  (Table  9) 
and  +0.17±0.07  for  the  corrected  data  (Table  10).  For  the  Tonga  earthquake  it  b  -H).07 
for  the  uncorrected,  -H).06  for  the  corrected  data.  As  noted  above,  the  Soviet  stations  are 
near  a  maximum  on  the  radiation  pattern  for  the  Mammoth  Lakes  earthquakes,  and  near 
the  null  axb  for  the  Tonga  event. 

Several  publbhed  works  treat  magnitude  reriduab  for  stations  of  the  Soviet  network, 
and  provide  the  basb  for  a  comparison  with  our  results  (Table  12).  First,  in  a  study  of 
magnitude  and  global  network  detection  capability  Ringdal  (1985)  recomputed  for 
about  70,000  earthquakes,  using  A  and  T  values  given  in  the  ISC  BvUetin  and  a 
maximum-likelihood  estimation  technique  (Ringdal,  1976).  In  another  study,  North 
(1976)  calculated  mean  station  magnitude  bias  from  values  given  in  the  ISC  Bulletin 
for  nearly  40,000  events  from  1964  to  1973.  The  biases  were  computed  for  the  "best"  (in 
terms  of  events  reported)  72  stations  with  respect  to  the  mean  magnitude  of  observations 
reported  by  thb  set  of  stations,  with  the  requirement  that  an  event  be  reported  by  more 


than  15  of  these  stations  before  a  bias  was  calculated. 

In  work  by  Soviet  authors  Vanek  et  »L  (1978,  1980)  determined  magnitude  correc- 
tioru  (Am^)  for  32  reference  stations  of  the  Unified  System  of  Seismic  Observations 
(ESSN)  of  the  Soviet  Union,  following  a  recommendation  of  the  KAPG  Conference  at 
Prague  in  1972  to  create  a  uniform  system  for  determining  magnitude  for  the  Eurasian 
continent.  The  station  corrections  were  determined  separately  for  phases  PV,  PVg  and 
PH  —  respectively  the  P-wave  recorded  on  vertical  mid-band,  vertical  short-period  and 
horisontal  mid-band  seismographs.  Calculations  were  also  made  separately  for  five  dif¬ 
ferent  source  regions  around  the  USSR  —  Alaska,  Japan,  the  Phillipines,  Asia  and  the 
Mediterranean.  The  number  of  earthquakes  in  each  source  region  was  not  given  in  the 
1980  paper,  but  in  the  mrUer  work  it  ranged  from  35  events  for  Asia  and  the  Mediter¬ 
ranean  to  93  events  for  Japan.  The  corrections  were  with  respect  to  one  of  the  reference 
stations,  OBN,  selected  at  least  in  part  for  its  bias  toward  large  values.  Based  on  a 
comparison  of  Am^  values  for  the  various  source  regions  the  reference  stations  are 
grouped  according  to  the  number  of  corrections  needed  for  the  different  source  regions. 
Thus,  a  station  for  which  Am^  was  the  same  for  all  five  source  regions  was  classed  as  a 
reference  station  of  the  /  Kind,  a  station  with  the  same  correction  for  four  of  the  five 
source  redone  would  be  a  station  of  the  //  Kind,  etc.  Of  interest  to  our  study,  station 
SEM  has  the  same  PVg  correction  for  all  of  the  source  regions,  Am^  »  +0.32,  making  it  a 
station  of  the  /  Kind  and  indicating  that  it  has  a  Sm^  bin*  of  -0.32  relative  to  station 
OBN.  Standard  deviations  are  not  given  by  Vanek  et  ni  (1980)  but  in  the  earlier  work 
they  average  ±0.05  unit.  For  comparison  with  magnitude  bias  given  by  other  authors, 
we  reversed  the  sign  of  Am^  and  increased  all  of  the  resulting  values  by  0.39  —  Ringdal’s 
(1985)  bias  for  station  OBN,  which  Vanek  et  cL  use  as  a  base  station. 

With  a  couple  of  exceptions  the  station  residuals  listed  in  Table  12  from  the  work  of 
Ringdal  (1985),  North  (1976)  and  Vanek  et  nL  (1978,  1980)  agree  to  within  a  few  hun¬ 
dredths  of  a  unit  of  m^.  Our  station  residuals  are  based  on  a  very  limited  data  set  and 
show  more  scatter  when  compared  to  those  of  the  other  authors.  However,  they  are  in 
general  agreement  with  the  published  values,  and  two  of  the  stations  —  SEM  and  YAK 
have  values  of  that  are  in  excellent  agreement  with  those  of  Vanek  et  oL  and  Ringdal, 
respectively. 

As  a  final  step  we  computed  the  bias  between  a  granite  site  at  the  Nevada  Test 
Site  and  the  Soviet  stations  listed  in  Table  12.  In  this  calculation  we  used  A  and  T  meas¬ 
urements  listed  by  Der  et  of.  (1978)  for  83  seismic  events  for  the  period  September  1976  to 
March  1977,  recorded  by  a  digital  seismic  system  (SDCS)  at  station  OB2-NV  (Climax 
stock  at  north  end  of  Yucca  Valley  on  NTS).  In  the  Der  et  oL  report,  the  A  values  are 
peak-to-peak  amplitudes  in  nanometers,  and  magnitudes  are  computed  using  distance 
corrections  of  Veith  and  Clawson  (1972).  For  consistency  with  the  work  by  Ringdal 
(1985)  and  North  (1976),  magnitudes  were  recomputed  using  the  Gutenberg  and  Richter 
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1  "  Giren  by  authon  comctiotu  rdatira  to  baaa  atatioD 

OBN.  Sign  rerenod  and  all  raluoa  incrcaiad  by  0.39  for  comparison  with 
Ringdal  (1925). 

2  "  Mean  of  Sm^  for  Tonga  eaithquako  phu  aTctago 

Stn^  for  four  Mammoth  Lako*  craiU  from  Table  9.  Uncoirectcd  data. 

3  "  Mean  of  Stn^  for  Tonga  eaithqnako  phu  areraga 

Stn^  for  four  Mammoth  Lakoa  ormts  from  Table  10.  Corrected  data. 

4  -  ATeiage  residual  (-0.10±0.35)  for  OB2-NV  with  respect  to  ISC 
BuUrtin  minus  aTeiage  station  residual. 

(1956)  corrections  and  dividing  the  peak-to-pesk  amplitudes  by  two  to  obtain  2ero-to- 
peak  amplitudes.  Residuals  for  the  83  events  were  calculated  as 
“  fn^(OB2-NV)-m^(ISC),  and  these  were  averaged  to  obtain  an  average  Sm.  of 


•0.09±0.39.  Five  values  that  fell  outside  the  2tr  limits  (—0.87  ^  ^  0.69}  were  dropped 

and  the  average  recalculated  to  obtmn  *0.10^0.35.  This  number  represents  the  aver^ 
age  bias  of  the  OB2-NV  site  with  respect  to  network-averaged  values  listed  in  the  ISC 
Bulletin.  To  determine  the  bias  of  OB2-NV  with  respect  to  the  Soviet  stations,  the  aver¬ 
age  residuals  in  Table  10  for  those  stations  were  subtracted  from  -0.10.  The  resulting  bias 
values  are  listed  in  the  right-hand  column  of  Table  12. 

Of  particular  interest  to  questions  of  yield  verification,  the  bias  of  the  NTS  gran¬ 
ite  site  with  respect  to  station  SEM  at  Senupalatinsk  is  -0.20,  with  a  range  of  -0.17  to 
-0.22.  The  smaller  of  these  figures  (-0.17)  is  based  on  the  study  by  Vanek  et  oL  (1980), 
and  the  larger  (-0.20,  -0.22)  are  from  our  measurements  of  Soviet  records  of  the  Mam¬ 
moth  Lakes  and  Tonga  earthquakes.  The  reader  should  be  reminded  that  the  Semipala- 
tinsk  station  is  about  100  km  from  the  Elast  Kasakh  test  site,  and  that,  according  to 
Peyve  and  Mbssakovsky  (1982),  crustal  rocks  under  the  test  site  have  been  subjected  to 
greater  folding,  faulting  and  intrusion  than  those  under  the  s^mic  station. 

It  should  also  be  noted  that  the  bias  values  in  Table  12  represent  only  the  bias  due 
to  attenuation  in  the  upper  mantle  and  crust  under  the  respective  seumic  stations;  they 
do  not  include  other  effects  such  as  differences  in  coupling  for  explosions  at  the  two  test 
sites,  effects  due  to  tectonic  release,  or  those  related  to  focusing  and  defocusing  of  sebmic 
waves  in  the  vicinity  of  a  given  explosion.  It  is  also  possible  that  magnitude  determina¬ 
tions  in  this  study  were  biased  because  of  the  different  passbands  of  the  instruments  used 
at  different  stations.  In  an  early  study  (SIPRI,  1968)  Whitham  reported  that 
intermediate-band  instruments  in  Canada  gave  magnitudes  that  averaged  about  0.3  m^ 
larger  than  short-period  instruments;  in  the  same  report  Karnik  reported  that  magnitudes 
based  on  broadband  Kimos  recordings  were  0.5  greater  than  for  a  Western  narrow- 
band  instrument,  but  the  difference  was  only  0.2-0.3  m^  if  a  short-period  SVK-M  instru¬ 
ment  was  used.  In  an  extension  of  this  study  we  plan  to  convolve  the  ground  motion 
illustrated  in  Figures  13-17  with  the  SDCS  instrument  response,  for  a  better  comparison 
with  the  magnitudes  determined  from  data  collected  at  NTS. 
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APPENDIX  A 


Desciipticm  of  the  SKM-S  System  and  Determination  of 
Response  Characteristics  for  the  Data  Used  in  this  Study 

According  to  Aranovich  et  mL  (1974),  "systems  of  [the  SKM-3]  type  are  used  pri¬ 
marily  for  the  recording  of  body  waves,  generated  by  earthquakes  of  medium  intensity 
(3  <  Af  <  5-6)  at  relatively  large  epicentral  distances  (greater  than  1,000  km).  For  this 
it  is  necessary  to  have  magnification  of  the  ordn  of  30,000-100,000  for  periods  of  ground 
motion  in  the  range  0.S-1.S  seconds.  Further  increasing  the  magnification  or  broadening 
the  frequency  response  to  longer  period  is  impoarible  without  special  filtering,  because  of 
microeuams  of  the  first  type  and  other  seismic  noise.  The  high  gains  mentioned  above 
can  be  achieved  only  at  stations  located  in  especially  favorable  surface  conditions  and  suf¬ 
ficiently  distant  firom  sources  of  interference.  Instrument  constants  and  gain  level  are  a 
function  of  the  noise  level  at  each  rite.  To  reduce  the  influence  of  local  arismic  noise  four 
standard  types  of  amplitude  response  are  used.  The  instrument  constants 

m,  p,  f,  s  and  U _ (where  U _ is  the  maximum  value  of  the  fluency  response) 

corresponding  to  the  four  types  of  response  are  given  in  Table  Al.  The  corresponding 
amplitude  response  curves  arc  shown  on  Figure  Al." 

The  standard  SKM-3  sosmograph  system  connate  of  two  SGKM-3  (horixontal)  and 
one  SVKM-3  (vertical)  seiamometers,  three  GK-VIIM  galvanometers  and  a  PS-3M  drum 
recorder.  The  latter  registers  the  light  beams  from  the  three  galvanometers  on  a  29-cm 
wide  by  90-cm  long  photographic  recording.  Time  marks  from  a  chronometer  are  printed 
once  per  minute.  Drum  speed  was  60  mm/nunute  for  all  of  the  SKM-3  stations  used  in 
this  study  except  KZL,  which  had  drum  speed  of  120  mm/minute. 

In  general,  for  a  seismometer-galvanometer  system,  the  relative  amplitude  response 
U  at  period  T  is  given  by  (Archangel’skii  et  at,  1974) 
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where  T  and  D  are  the  period  and  damping  factor,  respectively,  subscripts  a  and  g  refer  to 

2 

the  seismometer  and  galvanometer,  and  o  b  the  coupling  coefiicient. 

The  total  gun  of  the  system  is  given  by 
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where  A  is  the  optical  lever,  L  is  the  length  of  the  pendulum,  and  are  the  moments 
of  inertia  of  the  seismometer  and  galvanometer,  and  the  other  parameters  are  as  given 
above. 


The  phase  shift  at  period  T  is  given  by 


-1  +  p2^  -  *r* 
tan  Tf  = - - 

mT-q-f 


(4) 


Depending  on  the  period  T  the  phase  shift  may  vary  from - (for  T  =  0)  to  —  (for  T 

2  2 

-oo). 

Table  A2  lists  the  seismometer  and  galvanometer  constants  used  to  C2dculate 
response  curves  (Figure  6,  above)  for  the  instruments  used  in  this  study,  from  equations 
(1)  and  (2).  For  stations  BOD,  IRK  and  SEM  the  constants  were  taken  from  Shishkevish 
(1974)  and  the  response  curves  were  checked  to  insure  that  the  range  of  was  the  s£une 
as  indicated  on  the  recordings  or  in  the  Soviet  bulletins.  For  station  ELT  the  constants 
given  by  Shishkevish  were  modified  to  give  the  appropriate  upper  limit  of  I’m.  and  the 
same  constants  were  used  for  UST.  For  station  YAK  Shishkevish  lists  instrument  param¬ 
eters  only  for  the  horizontal-component  seismographs,  and  these  give  a  different  range  of 
than  that  on  the  records.  For  KZL  and  SEI  Shishkevish  does  not  list  instrument  con¬ 
stants.  As  a  result,  for  KZL,  SEI  and  YAK  various  combinations  of  instrument  constants 
were  tried,  based  primarily  on  tables  of  standard  setups  given  by  Shishkevish  and  by 
Aranovich  et  af.,  until  response  curves  matched  the  range  of  marked  on  the  records. 


Tsblc  AS.  bstrument  Parameters  for  Data  in  this  Study 
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Figure  Al.  Standard  amplitude  response  curves  for 
SKM-3  system  (Aranovich,  et  of.,  1974). 
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